DESY-03-059 
April 2003 



Measurement of deeply virtual Compton 
scattering at HERA 



ZEUS Collaboration 



Abstract 

The cross section for deeply virtual Compton scattering in the reaction ep —>■ ejp 
has been measured with the ZEUS detector at HERA using integrated luminosi- 
ties of 95.0 pb _1 of e + p and 16.7 pb _1 of e~p collisions. Differential cross sec- 
tions are presented as a function of the exchanged-photon virtuality, Q 2 , and the 
centre-of-mass energy, W, of the 7*p system in the region 5 < Q 2 < 100 GeV 2 
and 40 < W < 140 GeV. The measured cross sections rise steeply with increas- 
ing W. The measurements are compared to QCD-based calculations. 
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1 Introduction 



This paper reports cross-section measurements for the exclusive production of a real 
photon in diffractive ep interactions, ep —>■ eyp, as shown in Fig. This exclusive 
process, known as deeply virtual Compton scattering (DVCS) [HEIEHI] , is calculable in 
perturbative QCD (pQCD), when the virtuality, Q 2 , of the exchanged photon is large. 

The DVCS reaction can be regarded as the elastic scattering of the virtual photon off 
the proton via a colourless exchange. The pQCD calculations assume that the exchange 
involves two partons, having different longitudinal and transverse momenta, in a colourless 
configuration. These unequal momenta are a consequence of the mass difference between 
the incoming virtual photon and the outgoing real photon. The DVCS cross section 
depends, therefore, on the generalised parton distributions (GPD) [SHUIM!] , which carry 
information about the wave function of the proton. The cross section at sufficiently large 
Q 2 is expected to rise steeply with increasing W, the centre-of-mass energy of the virtual 
photon-proton system, due to the fast rise of the parton densities in the proton towards 
smaller x values, where x is the Bjorken scaling variable. 

The initial and final states of the DVCS process are identical to those of the purely 
electromagnetic Bethe-Heitler (BH) process (Figs. [B>,c). The interference between these 
two processes in principle provides information about the real and imaginary parts of the 
QCD scattering amplitude 00110] • However, it is expected to be small in the kinematic 
region studied in this paper jHJE]- 

The simplicity of the final state and the absence of complications due to hadronisation 
mean that the QCD predictions are more reliable than for many other exclusive final 
states. This reaction is one of the theoretically best-understood exclusive QCD processes 
in ep collisions. The first measurements of the DVCS process at high W fBEj arid its 
beam-spin asymmetry in polarised ep scattering at low W |13| IT4*] have recently become 
available. 

In the analysis presented here, the dependence of the DVCS cross section on W and Q 2 
is studied in the kinematic range 5 < Q 2 < 100 GeV 2 and 40 < W < 140 GeV. The 
measurements are integrated over t, the square of the four-momentum transfer at the 
proton vertex. The e + p cross sections are based on a ten-fold increase in statistics over a 
previous HERA result [12], permitting a study of the W dependence of the cross section 
as well as a significant extension of the Q 2 range probed. This paper also reports the first 
measurement of the e~p cross sections. 
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2 Experimental set-up 



The data were collected by the ZEUS detector at HERA during the 96-00 running periods. 
In 96-97, HERA collided 27.5 GeV positrons with 820 GeV protons. In 98-00, the proton 
energy was increased to 920 GeV and both positrons and electrons were collided. The 
measurements for e + p (e~p) interactions 1 are based on an integrated luminosity of 95 pb _1 
(^pb- 1 ). 

A detailed description of the ZEUS detector can be found elsewhere jT^]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) ^B], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organised in nine superlayers covering 
the polar-angle 2 region 15° < 6 < 164° (2.02 > rj > —1.96). The transverse-momentum 
resolution for full-length tracks is cr{p T )/p T = 0.0058pr © 0.0065 © 0.0014/p r , with p T in 
GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) consists of three parts: 
the forward (FCAL, 1.1 < 77 < 3.8), the barrel (BCAL, -0.7 < rj < 1.1) and the rear 
(RCAL, —3.4 < 7] < —0.7) calorimeters. Each part is subdivided transversely into towers 
and longitudinally into one electromagnetic section (EMC) and either one (RCAL) or two 
(BCAL and FCAL) hadronic sections (HAC). The smallest subdivision of the calorimeter 
is called a cell. In the EMC section, the towers are divided transversely into either four 
(FCAL and BCAL) or two (RCAL) cells. The CAL energy resolutions, as measured under 
test-beam conditions, are a(E)/E = 0.18 /y/E for electrons and a(E)/E = 0.35/ \/E for 
hadrons, with E in GeV. 

Presampler detectors ^H] are mounted in front of the CAL. They consist of scintillator 
tiles matching the calorimeter towers and measure signals from particle showers created 
by interactions in the material lying between the interaction point and the calorimeter. 
In this analysis, only the information from the presampler in front of the RCAL was used 
to correct the energy of the final-state particles. 

The forward plug calorimeter (FPC) ^H] is a lead-scintillator sandwich calorimeter with 
readout via wavelength-shifter fibres. It was installed in 1998 in the 20 x 20 cm 2 beam 

1 Hereafter, both e + and e~ are referred to as electrons, unless explicitly stated otherwise. 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as rj = — In (tan |), where the polar angle, 6, is measured with respect to the proton beam 
direction. 
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hole of the FCAL and has a small hole of radius 3.15 cm in the centre to accommodate the 
beam pipe. It extends the pseudorapidity coverage of the forward calorimeter from rj < 4.0 
to 77 < 5.0. The FPC information was used to remove low-mass proton-dissociative events 
from the analysis. 

The hadron-electron separator (HES) (201 i s installed in the RCAL and FCAL. It consists 
of 3 x 3 cm 2 silicon diodes placed at a longitudinal depth of three radiation lengths, which 
corresponds to the approximate position of the maximum of the electromagnetic shower 
in the CAL. The separation between electrons and hadrons is based on the fact that the 
hadronic interaction length is 20 times larger than the electromagnetic radiation length. 
In this analysis, the fine segmentation of the RHES was used to improve the position 
resolution for both scattered electrons and photons. 

The small-angle rear tracking detector (SRTD) |2^ is attached to the front face of the 
RCAL {Z = —148 cm). The SRTD consists of two planes of scintillator strips read out 
via optical fibres and photomultiplier tubes. It covers the region 68 x 68 cm 2 in X and 
Y with the exclusion of a 8 x 20 cm 2 hole at the centre for the beam pipe. The SRTD 
provides a transverse-position resolution of 3 mm and was used to measure the positions 
of photons and electrons scattered at small angles relative to the lepton beam direction. 

The proton-remnant tagger (PRT1) [22] consists of two layers of scintillation counters 
located at Z = 5.15 m, and covers the pseudorapidity range 4.3 < r] < 5.8. It was used, 
up to the end of the 1997 running period, to tag events in which the proton diffractively 
dissociated. 

The luminosity was determined from the rate of the bremsstrahlung process ep — > ejp, 
where the high-energy photon was measured with a lead-scintillator calorimeter [23] lo- 
cated at Z = —107 m. 

3 Event selection 

For the Q 2 range of this analysis, Q 2 > 5 GeV 2 , and small t, the signature of DVCS 
and BH events consists of a photon and a scattered electron with balanced transverse 
momenta. The scattered proton stays in the beam pipe and remains undetected. 

The events were selected online via a three-level trigger system [T5JI2I]- The trigger 
selected events with two isolated electromagnetic (EM) clusters in the EMC with energy 
greater than 2 GeV. The events were selected offline by requiring two EM clusters, the 
first in the RCAL with energy E\ > 15 GeV and the second, with polar angle 0.6 < 9 2 < 
2.75 rad (1.2 > 772 > -1.6), either in the RCAL, with energy E 2 > 3 GeV, or in the BCAL, 
with energy E 2 > 2.5 GeV. The angular range of the second cluster corresponds to the 
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region of high efficiency for reconstruction of a track in the CTD. if a track was found, 
it was required to match one of the EM clusters. Events with more than one track were 
rejected. To ensure full containment of the electromagnetic showers in the CAL, events 
in which one of the clusters was located within 3 cm of the beam hole were rejected. 

The selection 40 < E — pz < 70 GeV was imposed, where E is the total energy and pz 
the sum of EcosO over the whole CAL. This requirement rejects photoproduction events 
and events in which a hard photon is radiated from the incoming electron. 

After these cuts, when the two EM clusters are ordered in energy such that E\ > E 2 , the 
kinematics ensure r]i < 772- In the following, the two clusters will be denoted as EMI and 
EM2, respectively. 

For the 96-97 (98-00) running period, calorimeter cells not associated with the two elec- 
tromagnetic clusters were required to have energy less than: 150 (200) MeV in the FEMC 
and 200 (300) MeV in the FHAC; 200 (350) MeV in the BEMC and 250 (350) MeV in the 
BHAC; 150 (150) MeV in the REMC and 300 (300) MeV in the RHAC. These thresholds 
were set to be three standard deviations above the noise level of the CAL. Moreover, for 
the 98-00 data sample, the energy measured in the FPC was required to be less than 
1 GeV. These elasticity requirements reject most events in which the proton dissociates 
into a hadronic system, X. 

The events were subdivided into three samples, a first in which there was no track associ- 
ated with the EM2 cluster (7 sample), a second in which the track associated with EM2 
had the same charge as the beam electron (e sample) and finally a third in which the track 
associated with EM2 had the opposite charge to that of the beam electron (wrong-sign-e 
sample). These samples are interpreted as: 

• 7 sample: EM2, with no track pointing to it, is the photon candidate and EMI is 
the scattered-electron candidate. Both BH and DVCS processes contribute to this 
topology. The sample consisted of 3945 events. 

• e sample: EM2, with the right-charge track pointing to it, is the scattered-electron 
candidate and EMI is the photon candidate. This sample is dominated by the BH 
process. The contribution from DVCS is predicted to be negligible, due to the large 
Q 2 required for a large electron scattering angle. This sample contained 7059 events. 

• wrong-sign-e sample: EM2, with the wrong-charge-sign track pointing to it, may have 
originated from an e + e~ final state accompanying the scattered electron, where one 
of the right-sign electrons escaped detection. This background sample is due to non- 
resonant e + e~ production and to J ftp production and subsequent decay. Other sources 
are negligible, as will be discussed later. This sample consisted of 287 events. 

The wrong-sign-e sample was used to statistically subtract the background contributions 
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to the e sample in each kinematic bin. The background-subtracted e sample was then 
used to investigate the BH contribution to the 7 sample. 

For the purposes of this analysis, the values of Q 2 and W were determined for each 
event, independently of its topology, under the assumption that the EMI cluster is the 
scattered electron. This assumption is always valid for DVCS events for the Q 2 range 
considered here. The value of Q 2 was calculated using the electron method [22], while 
W was determined using the double-angle method j2S|. No explicit cut on t was applied 
in the event selection. Events for which 40 < W < 140 GeV and 5 < Q 2 < 100 GeV 2 were 
retained. 



4 Monte Carlo simulations 

The acceptance and the detector response were determined using Monte Carlo (MC) simu- 
lations. The detector was simulated in detail using a program based on GEANT 3.13 |26j . 
All of the simulated events were processed through the same reconstruction and analysis 
chain as the data. 

A MC generator, GenDVCS [2Zj based on a model by Frankfurt, Freund and Strikman 
(FFS) |3], was used to simulate the elastic DVCS process. In the FFS calculation, the 
DVCS cross section, integrated over the angle between the e and p scattering planes, is 
related to the inclusive structure-function F 2 through 

J3 ep^e 7 p 2 3 

where x ~ Q 2 / (Q 2 + W 2 ) is the Bjorken scaling variable, b is the exponential slope of the 
t dependence and y is the fraction of the electron energy transferred to the proton in its 
rest frame. The ratio R = ^smA(j*p — > r y*p)\t=o/'^ m -^(l*P ~* lP)\t=o accounts for the 
non-forward character of the DVCS process and is directly related to a ratio of the GPD 
to the parton distribution functions [21] and p is the ratio of the real to imaginary part 
of the amplitude, p = WteA^p — > r yp)\ t= o/^smA( , j*p — > 7p)|t=o- 

The value of R, calculated using the leading-order (LO) QCD evolution of the GPD, is 
about 0.55, with little dependence on x or Q 2 jlj. For simplicity, the R parameter in the 
MC generator was set to a constant value of R = 0.55. 

In GenDVCS, the ALLM97 [2E] parameterisation of the F2 structure function of the 
proton was used as input. In this empirical fit to 7*p total-cross-section data, the value 
of p was parameterised as p = |(0.176 + 0.033 In Q 2 ), where Q 2 is in GeV 2 [21!. 

In the FFS model, the t dependence is assumed to factorise, with the slope parameter, b, 
depending on both W and Q 2 . The value of b is expected to decrease with Q 2 and, even 
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at high Q 2 and at very small x, is expected to increase with W. While this dependence is 
important for the normalisation of the calculated DVCS cross section, it does not affect 
the acceptance corrections. In the MC simulation b was assumed to be constant and was 
set to 4.5 GeV" 2 . 

For proper treatment of radiative effects, the GenDVCS generator was interfaced to HER- 
ACLES 4.6 [20], which includes corrections for initial- and final-state photon emission from 
the electron line, as well as vertex and propagator corrections. 

The elastic and inelastic BH processes, ep — > e^p and ep — > e^X, and the QED di- 
lepton production, ep — > ee + e~p, were simulated using the GRAPE-Compton 3 [Slj and 
GRAPE-Dilepton generators, respectively. These two MC programs are based on the 
automatic system GRACE [22! for calculating Feynman diagrams. The GRAPE generator 
gave identical results to the Compton 2.0 jHS] generator for the elastic BH process. The 
GRAPE program was used because it simulates the hadronic final state for the inelastic 
BH process. 

Additional samples were generated using the diffractive RAPGAP and non-diffractive 
DJANGOH (33] generators in order to study possible backgrounds from low-multiplicity 
DIS events. A possible contribution from vector-meson electroproduction was simulated 
by the ZEUSVM [H] MC generator interfaced to HERACLES. 

5 DVCS-signal extraction 

In the kinematic region of this analysis the interference between the DVCS and BH 
amplitudes is very small when the cross section is integrated over the angle between 
the e and p scattering planes [HUH]- Thus the cross section for exclusive production of 
real photons may be treated as a simple sum over the contributions from the DVCS and 
electromagnetic BH processes. The latter can, therefore, be subtracted and the DVCS 
cross section determined. 

The BH process was studied using the e sample which, according to the MC predic- 
tions, consists almost solely of BH candidates. A background contribution to the BH 
events in this sample of about 4%, originating from deep inelastic exclusive e + e~ pro- 
duction, where one of three final-state leptons escapes detection, was estimated from the 
wrong-sign-e sample and MC simulations, and was statistically subtracted. According to 
the MC simulations, 75% of this background consists of non-resonant di-lepton production 
and 25% of exclusive J/ip production with subsequent e + e~ decay. The normalisation of 

3 Hereafter, the GRAPE-Compton generator is referred as GRAPE 
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the MC samples was determined from the wrong-sign-e sample. The diffractive electro-pro- 
duction of p, u and mesons, in which one of the decay charged particles was misidentified 
as an electron in the CAL, the other was undetected, and the electron scattered into the 
RCAL was taken to be the photon, was negligible. 

The expectation for the inelastic BH contribution to the e sample is subject to un- 
certainties coming from the dependence of the selection efficiency on the mass of the 
hadronic final-state system, X. This inelastic contribution was estimated from the data 
as (17.8 ± 1.2)%, for the 96-97 data sample, using the fraction of events tagged in the 
PRT1 and (10.5 ± 1.0)%, for the 98-00 data sample, using events with more than 1 GeV of 
energy in the FPC, obtained releasing the elasticity cut. The uncertainties are statistical. 
No attempt was made to quantify the systematic uncertainty since there is little sensi- 
tivity to the cross section. The difference in the measured fraction of proton-dissociative 
events is due to differences in the detector configuration between the two running periods. 

After the subtraction of the dilepton and J /if) backgrounds, for the 96-97 (99-00) e + p 
samples, the number of remaining BH events in the e sample was 2523 (3289), while the 
expected number from the GRAPE simulation was 2601 (3358). The absolute expectation 
of the GRAPE simulation reproduced the number of BH data events to within (3 ± 5)% 
for 96-97 and (2 ± 4)% for 99-00, where uncertainties include the statistical uncertainty 
as well as the uncertainties due to the trigger efficiency and the estimation of the inelastic 
BH contribution. 

For further analysis, the GRAPE MC sample was normalised to the BH events in the 
data. The comparison between the e sample and the sum of the GRAPE BH, dilepton 
and J/ifj MC samples for the 99-00 e + p running period is shown in Fig. |2K versus W and 
in Fig.Eb versus the difference in azimuthal angles of the EMI and the EM2, A0 12 . Good 
agreement is observed. 

The properties of the 7 sample were then compared to the expectations of the normalised 
GRAPE MC sample. As an example, the comparison of the W and A0i 2 distributions 
for the 99-00 e + p running period is shown in Figs. and |21i, respectively. An excess of 
events over the expectations of the GRAPE simulation is observed. Moreover, the data 
distributions differ from those expected for the BH process. The W distribution of the 
BH sample peaks at large W, while that of the data is more evenly distributed. The A</> 12 
distribution of the BH sample is also narrower than that of the data. 

The W and the A0i 2 distributions of the data, after subtracting the BH contribution 
using the renormalised GRAPE MC sample, are shown in Figs. and I2F, respectively. 
This sample includes events in which the proton dissociated into a hadronic final state 
with low mass. 

The previous ZEUS measurements of elastic vector meson production j^TJEH] support, 
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within relatively large uncertainties, the assumption that the fraction of proton dissocia- 
tive events, / p _diss, in diffractive interaction is process independent. Therefore, in this 
analysis, the values of /p-diss determined from the measurements of the diffractive J/psi 
photoproduction [HB] are used 

/p-diss = 22.0 ± 2.0(stat.) ± 2.0(syst.)% for the 96-97 data; 
/ p _diss = 17.5 ± 1.3(stat.)t|£(syst.)% for the 98-00 data. 

The above fractions are consistent, within large uncertainties, with those estimated using 
the events in the 7 sample either tagged by the PRT1 or the FPC, after subtracting the 
inelastic BH contributions. 

Other possible sources of contamination were investigated. Due to the relatively high Q 2 
of the present data set, the contamination from production of light vector mesons, such as 
uj or 0, decaying through channels containing photons in the final state is below 1% and 
was neglected. A possible contribution to the 7 sample from low-multiplicity processes 
such as ep — > eir p, ep — > e7r°7r°p, and ep — > eii rip, where the tt° or 77 fakes a photon 
signal, was also investigated. The number of candidate events found in the RAPGAP and 
D JANGOH samples was reweighted to reflect the cross sections obtained by extrapolating 
low- IV measurements pMl4Ti] . Their contribution to the DVCS sample is negligible. 

The data are compared to the absolute expectations of GenDVCS in Figs. Eh andEF- The 
best agreement in normalisation between the data and the M C simulation is achieved when 
the normalisation of the latter is decreased by 10%. This was obtained by increasing the 
value of b from 4.5 to 4.9 GeV -2 (see Section HJ). Overall, good agreement between the 
data and the MC simulation is found, demonstrating that the excess of photon candidates 
over the expectation of BH is due to DVCS. 



6 Cross-section determination 

The 7*^ cross section for the DVCS process as a function of W and Q 2 was evaluated 
using the expression 

a( 7 > - 7P)(Wi, Qi) = (ATbS = " /P " di8s) • * FFS (7> - 7P)(Wi, <?■), 

where N° is the total number of data events in the 7 sample in bin % in W and Q 2 , 
iVf 11 denotes the number of BH events in the 7 sample in that bin, determined from the 
renormalised GRAPE sample, and is the number of events expected in the 7 sample 
from GenDVCS for the luminosity of the data. The factor / p _diss is the fraction of the 
proton-dissociative DVCS events in the data, a FF (*y*p — > 7p) is the j*p cross section 
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computed according to the FFS expression, and Wi and Q 2 are the values of W and Q 2 
where the cross section is evaluated. 

The 7*p cross sections have been computed in the ranges 5 < Q 2 < 100 GeV 2 and 40 < 
W < 140 GeV, separately for the 96-97, 98-99 and 99-00 data periods and then combined 
for the positron samples (96-97 and 99-00). Tables [T] - El list the j*p —> jp cross-section 
values. 

7 Systematic uncertainties 

The systematic uncertainties of the measured cross sections were determined by changing 
the selection cuts or the analysis procedure in turn and repeating the extraction of the 
cross sections. The following systematic studies have been carried out: 

• all the selection cuts discussed in Section |3] were shifted according to the resolutions 
of the corresponding variables. The most significant contributions came from varying 
the lower Q 2 cut. The average change in the cross section due to this cut was ±2%. 
The largest change in the cross section, ±10%, was found in the highest-l^ bin, while 
it was ±4% in the lowest-Q 2 bin; 

• the elasticity cut was changed by ±30 MeV in the EMC and ±50 MeV in the HAC. 
The average change in the cross section was ±2% in all bins of Q 2 and W, while the 
largest change in the cross section, observed when the cut was lowered, was —4% in 
the lowest-W bins and —4% in the lowest-Q 2 bin; 

• the trigger efficiency was varied within its statistical uncertainty. This resulted in 
average changes of the cross section of about ±2%. The biggest variation of the cross 
section of ±3% was observed in the lowest-Q 2 bin and in the two highest-VT bins; 

• the electromagnetic energy scale was varied within its uncertainty of 1.5% for the EM2 
(low energy) and of 1% for the EMI (high energy), resulting in a ±3% average change 
of the cross section in both Q 2 and W. The largest change was ±3% for the lowest-Q 2 
bins and ±5% for the highest-VT bin; 

• in GenDVCS, the Q 2 dependence was modified by introducing a Q 2 -dependent t slope 
using the formula b = 8(1 — 0.15 ln(Q 2 /2)) GeV -2 (see Section^)). The average change 
in the cross section was ±1%, with the largest variation of ±3% in the highest-Q 2 bin. 

The uncertainty on the proton-dissociative contribution, / p _diss, leads to an overall nor- 
malisation uncertainty of ±4.0% and ±3.5% for the e~p and e + p data, respectively. 

The systematic uncertainties typically are small compared to the statistical uncertainties. 
The individual systematic uncertainties, including that due to /p-dissj were added in 



9 



quadrature separately for the positive and negative deviations from the nominal cross- 
section values to obtain the total systematic uncertainties listed in Tables [D-El An overall 
normalisation uncertainty in the luminosity determination of ±1.8% and ±2.0% for the 
e~p and e + p data, respectively, was not included because it was small with respect to the 
above contributions. 

8 Results 

The W dependence of the DVCS cross section, ctdvcs = &{l*P — > 7P), for Q 2 = 9.6 GeV 2 
is shown in Fig. |3J separately for e + p and e~p interactions. Due to the limited statistics, 
the e~ sample is only shown in three W bins. There is agreement between the two samples. 

A fit of the form odvcs oc W 5 was performed separately for the positron and electron 
data. For the e + p data, the value 5 = 0.75 ± 0.15(stat.)^Q;Qg(syst.) is comparable to that 
determined for J/ip electroproduction [HB]- This steep rise in the cross section is a strong 
indication of the presence of a hard underlying process. The same fit to the e~p data 
yields 5 = 0.45 ± 0.36(stat.)+o;Qy(syst.), which is compatible with the e + p result. 

The positron sample has been further subdivided into three Q 2 ranges. The W dependence 
of Odvcs i n these three Q 2 bins is presented in Fig. 0] The results are compatible with 
no dependence of 5 on Q 2 although also with the increase with Q 2 observed in exclusive 
production of light vector mesons jHJEIj. 

The Q 2 dependence of odvcs; f° r W = 89 GeV, is shown in Fig. again separately for 
e + p and e~p interactions. There is no significant cross-section difference between the e + 
and e~ data, which is consistent with the assumption that the present measurement is 
insensitive to the interference term. 

A fit of the form Q~ 2n to the e + p data gives a value of n = 1.54±0.07(stat.) ±0.06(syst.). 
This value is lower than n ~ 2 which is characteristic of exclusive vector-meson produc- 
tion (23112]. The fit to the e~p data gives n = 1.69 ± 0.21(stat.)+°;^(syst.). 

9 Comparison with models 

In the presence of a hard scale (Q 2 » Aq CD ), the DVCS amplitude factorises into a 
hard-scattering coefficient, calculable in pQCD, and a soft part which can be absorbed 
in the GPD jH]. The kernels of the evolution equations for the GPD are known to next- 
to-leading order (NLO) and the GPD can thus be evaluated at all Q 2 given an 
input at some starting scale. Measurements of the DVCS cross section are an essential 
ingredient in modelling the input GPD |4*4*ll2"5] . 
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Freund, McDermott and Strikman (FMS) |2H] have made an attempt to model the GPD 
based on DVCS data (T^EUE] • A comparison of cxdvcs as a function of Q 2 for fixed 
W with the predictions based on the MRST parameterisation of the parton distribu- 
tion functions (PDF) |4^l4T)] is shown in Fig. Eb- Three FMS curves are shown. Two 
curves show the results of modelling based on LO (MRSTL) and NLO (MRSTM) parton 
distribution functions. The latter leads to predictions closer to the data. In this compar- 
ison, a fixed value of the t slope, b = 4.9 GeV -2 , was assumed. The third curve, shown 
in the figure, corresponds to predictions based on MRSTM, assuming a Q 2 -dependent 
b value. The best agreement between the data and the predictions is achieved using 
b = 8(1 — 0.15 ln(<5 2 /2)) GeV -2 , a parameterisation obtained by Freund, McDermott and 
Strikman j2H] from a fit to a preliminary version of the present data. Similar conclusions 
are reached when the CTEQ6 parameterisations [5711313132] are usec ^ ( no ^ shown). 

The data are also compared to the expectations of FFS (see Section |3|), again assuming 
b = 4.9 GeV -2 . For Q 2 > 20 GeV 2 , the e + p data lie significantly above the prediction. 

The DVCS cross section has also been calculated within colour-dipole models |oTIll51llo^l 
E21E31, which have been successful in describing both the inclusive and the diffractive 
DIS cross sections at high energy (HSlESlEZlEHlEniinni • The various dipole models differ 
in their formulation of the dipole cross section with the target proton. If s-channel 
helicity is conserved in DVCS, the virtual photon must be transversely polarised. As 
the wave function of the transversely-polarised photon can select large dipole sizes, whose 
interactions are predominantly soft, DVCS constitutes a good probe of the transition 
between perturbative and non-perturbative regimes of QCD. The Q 2 dependence of odvcs 
has been compared to the expectations of three calculations based on colour-dipole models, 
by Donnachie and Dosch (DD) [BU], Forshaw, Kerley and Shaw (FKS) jHTJEUISSl and 
McDermott, Frankfurt, Guzey and Strikman (MFGS) [B2l59j . The comparisons are shown 
in Fig. |H1 where the model expectations are represented by curves corresponding to a fixed 
value of b = 4 GeV -2 (upper) and b = 7 GeV -2 (lower), chosen for illustration. All three 
predictions give a reasonable representation of the data. The HI measurements ^2] are 
also shown, extrapolated to the W value of the ZEUS data using the W 5 dependence of the 
cross section measured in this analysis for the e + p data. The HI data lie systematically 
below the ZEUS data. 



10 Conclusions 

The DVCS cross section u{^*p — > jp) has been measured at HERA in the kinematic 
range 5 < Q 2 < 100 GeV 2 and 40 < W < 140 GeV. No significant difference between the 
e + p and e~p interactions was observed. The data have been compared to calculations 
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based on generalised parton distributions and on the colour-dipole model. Generally, 
good agreement with the data is observed. 

The Q 2 dependence of the DVCS cross section follows approximately a Q~ 3 behaviour. 
The precision of the data allows an accurate determination of the W distribution for the 
first time. The cross section rises steeply with W, indicative of a hard underlying process, 
where the rise reflects the increase of parton distributions with decreasing Bjorken x. 

These measurements demonstrate the potential of DVCS data to constrain the structure 
of the proton and quark-gluon dynamics at low x. 
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Q 2 range (GeV 2 ) 


Q 2 (GeV 2 ) 


e + p 


a 7>-7P ( n b) 
e~p 


5 - 10 


7.5 


5.42±0.33±g;§| 


5.63±0. 77^;| 


10 - 15 


12.5 


2.64±0.22lg;i3 


2.20±0.52±°;}| 


15 - 25 
25 - 40 
40 - 70 
70 - 100 


20 
32.5 

55 
85 


1.23±0.14i°$ 
0.59±0.12±g;g| 

0.20±0.08+°:o2 
0.16±0.09^o3 


0.96±0.3lig;£jj 
0.61±0.28t{j;gjj 



Table 1: Values of the cross sections for the 7*p — > 7p DVCS process as a 
function of Q 2 for the e + p and e~p data. Values are quoted at the centre of each 
Q 2 bin and for the average W value of the whole sample, W = 89 GeV , obtained 
from GenDVCS. The first uncertainty is statistical and the second systematic. The 
systematic uncertainty due to the luminosity determination is not included. 





W range (GeV) 


W (GeV) 


ct 7*p->7p (nb) 


W range (GeV) 


W (GeV) 


ct 7*p-7p (nb) 


e + p 


e + p 


e + p 


e~p 


e~p 


e~p 


40 - 50 


45 


2.19±0.24i°;^ 








50 - 60 
60 - 70 
70 - 80 


55 
65 
75 


2.96±0.28±g;i| 
3.62±0.36lg;^ 
3.88±0.42i^ 


40 - 73 


56.7 


2.94±0.39±g;i| 


80 - 90 


85 


3.59±0.45j££! 








90 - 100 


95 


3.29±0.55±g;|5 


73 - 107 


90 


4.06±0.69±g:ii 


100 - 110 


105 


6.24±0.77^jg 








110 - 120 


115 


4.86±0.76lg:44 








120 - 130 
130 - 140 


125 
135 


4.69±0.82+°;| 
5.55±0.99lg;|5 


107 - 140 


123.3 


3.8±l.l^i 



Table 2: Values of the cross sections for the 7*p — > 7p DVCS process as a 
function of W for the e + p and e~p data. Values are quoted at the centre of each 
W bin and for the average Q 2 value of the whole sample, Q 2 = 9.6 GeV 2 , obtained 
from GenDVCS. The first uncertainty is statistical and the second systematic. The 
systematic uncertainty due to the luminosity determination is not included. 
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(J 7*P- > 7P 


W range (GeV) 


W (GeV) 


a 7>-7P ( n b) 
5 < Q 2 < 8 GeV 2 
Q 2 = 6.2 GeV 2 


a 7>-7P ( n b) 
8 < Q 2 < 13 GeV 2 
Q 2 = 9.9 GeV 2 


a^P^P (nb) 
13 < Q 2 < 30 GeV 2 
Q 2 = 18.0 GeV 2 


40 - 65 
65 - 90 
90 - 115 
115 - 140 


52.5 
77.5 
102.5 
127.5 


5.63d=0.58±g;fg 
6.57±0.91±g;|J 
9.5±1.5±S:1 
7.6=hl.6±S:i 


2.52±0.26±g;?i 
3.12±0.39l°; 2 ^ 
3.94±0.61±g;^ 
5.83± 0.89±g;t| 


0.99±0.13j£So 

i.34±o.i7±8:8| 

1.91±0.30±g;i| 
1.64±0.47±g;il 



Table 3: Values of the cross sections for the 7*p — > 7p DVCS process as a 
function of W for the e + p data in three Q 2 ranges. Values are quoted at the centre 
of each W bin and for the average Q 2 values obtained from GenDVCS. The first 
uncertainty is statistical and the second systematic. The systematic uncertainty 
due to the luminosity determination is not included. 



17 




18 



0) 



® 10 3 



LU 



10 



10 



10 



10 



10 



10 



10 



10 



ZEUS 




ZEUS 99-00 e + p 
(e-sample) 
dilepton+JA]/ 



(a) 



■ ZEUS 99-00 e + p (c) 
(y-sample) 



(e) 

o ZEUS 99-00 e + p 

(y-sample BH subtracted) 



BH (GRAPE) 
BH+dilepton+JA|/ 



(b) 



BH+DVCS 



DVCS (GenDVCS) (f) 



60 80 100 120 0.95ti n 1.05ti 

W (GeV) A<|> 12 (rad) 



Figure 2: Comparison of the W and the A0 12 distributions for the data from 
99-00 e + p running period and corresponding MC samples as described in the figure; 
(a) and (b) e sample compared to the sum of the BH, dilepton and J ftp MC simu- 
lations, with the sum of the latter two also shown separately; (c) and (d) 7 sample 
compared to the expectation of the BH MC simulation, normalised to the data in 
the e sample after background subtraction; (e) and (f)'y sample after subtracting 
the BH expectation compared to GenDVCS, where the absolute normalisation of the 
latter is increased by 17.5% to account for the proton-dissociative component in the 
DVCS data. 
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Figure 3: The DVCS cross section, er(7*p — > 7p) ; as a function of W for an 
average Q 2 = 9.6 GeV 2 , separately for e + p data (dots) and e~p data (triangles). 
The solid line is the result of a fit of the form er DVCS ^ W 5 to the positron data. 
The error bars denote the statistical uncertainty (inner) and the quadratic sum of 
the statistical and the systematic uncertainties (outer). 
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Figure 4: The DVCS cross section, <r(7*p — > 7p), as a function of W for three 
Q 2 values for e + p data as denoted in the figure. The corresponding ranges in Q 2 
are listed in Tabled The solid line is the result of a fit of the form cr DV cs W 5 . 
The values of 5 and their statistical uncertainties are given in the figure. The last 
data point for Q 2 = 9.9 GeV 2 is displaced horizontally for ease of visibility. The 
error bars denote the statistical uncertainty (inner) and the quadratic sum of the 
statistical and the systematic uncertainties (outer). 
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Figure 5: (a) The DVCS cross section, a(j*p — > jp), as a function of Q 2 for 
W = 89 GeV, separately for e + p data (dots) and e~p data (triangles). The solid 
line is the result of a fit of the form odvcs ^ Q~ 2n to the positron data. The 
e~p data points are displaced horizontally for ease of visibility; (b) <r(7*p — > jp) 
as a function of Q 2 compared to the GPD-based theoretical predictions of FFS and 
FMS, where MRSTL(M) indicates the LO (NLO) parameterisation of PDF. The 
MRSTM expectations are also shown for the Q 2 -dependent b values described in 
the text. The error bars denote the statistical uncertainty (inner) and the quadratic 
sum of the statistical and the systematic uncertainties (outer). 
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Figure 6: The DVCS cross section, a(j*p — > jp), as a function of Q 2 , as also 
shown in Fig.\^ for the ZEUS (dots) and HI (squares) data. The data are compared 
to the theoretical predictions of the DD, FKS and MFGS models of colour-dipole 
interactions. The curves correspond to fixed b values, b = 4 GeV~ 2 (upper) and 
b = 7 GeV~ 2 (lower). 
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